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ABSTRACT 

We use the Eighth Data Release of Sloan Digital Sky Survey (SDSS DR8) catalog with its additional 
sky coverage of the southe rn Galactic hemi sphere, to measure the extent and study the nature of the 
Virgo Overdensity (VOD: Uuric et al1l2008f ). The data show that the VOD extends over no less than 
2000 deg^, with its true extent likely closer to 3000 deg^. We test whether the VOD can be attributed 
to a tilt in the stellar halo ellipsoid with respect to the plane of the Galactic disk and find that the 
observed symmetry of the north-south Galactic hemisphere star counts excludes this possibility. We 
argue that the Virgo Overdensity, in spite of its wide area and cloud-like appearance, is still best 
explained by a minor merger. Its appearance and position is qualitatively similar to a near peri- 
galacticon merger event and, assuming that the VOD a nd the Virgo Stellar S t ream share the same 
progenitor, consistent with the VSS orbit determined bv iGasetti-Dinescu et al.l ()2009[ ). 
Subject headings: Galaxy: structure — Galaxy: halo — Galaxy: formation 



1. INTRODUCTION 

A growing number of observational campaigns over the 
last decade showed the richness of substructure in the 
Milky Way halo. The most vivid discoveries of ongoing 
accretion onto the Galaxy include the Sagi ttarius dwarf 
and its tidal streams (Ivezic et al. 2000 t lYannv et all 
120001 iVivas et al.l [20011 : IMaiewski et al. '2003', amoi^ 
others ), the Monoceros stellar stream (.Ncwbcr g et al.l 
[2OOI iRocha-Pinto et al.l [200l. as weh as numerous 
smaller streams and dwarf galaxies abundant in the so- 
called "field of streams" (|I3elokurov et al.l [20061 ) . Ev- 
idence of ongoing mergers has also be en detected in 
the halo of M31 (gbata et al. 2001; Ferg uson et al.ll2002l: 
'Morrison et al.' 2003; Zucker et al.l 120041 : iKalirai et al.l 
2006. : Fardal et al., ,2007[), providing further support for 
the importance of hierarchical merging in galaxy forma- 
tion. 

Even if one is ultimately interested in the processes of 
galaxy formation in general, the study and understand- 
ing of these processes in the Milky Way remains advan- 
tageous in many aspects. For example, the proximity of 
stars and substructures in the Galaxy allows data col- 
lection at high resolution, which can be directly tested 
against predictions from simulations. 

In this wo rk, we concentra te on the Virgo Overden- 
sity (VOD; lJuric et al.l 120081 ) . a cloud-like overdensity 
of stars spanning distances between 10 — 20 kpc. Parts 
of the VOD were initially identified as an overden- 
sity of RR Lyrae stars confined to a small region in 
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the Q uasar Equatorial Survey Team survey (jVivas et al.l 
l2001f ). followed by an identification as an overdensity 
of main sequence tu rnoff stars in the same direction 
(jNewberg et al.l r2002| ) . However, it was not until the 
advent of wide- area surveys such as the Sloan Digital 
Sky Survey (SDSS; lYork et al.l [200(TI) that the size o f 
the VOD feature was fully recognized puric et al.ll2008() . 
In particular, u sing a large photo metric sample of main 
sequence stars, I Juric et "all (|2008[ ) have shown that the 
VOD extends over at least 1000 deg^, spans heliocentric 
distances between 6 < D < 20 kpc, and also shows as 
an overdensity of M giants in the Tw o Micron All Sky 
burvey (2MASS; iSkrutskie et~alll2006( ). Based on their 
u — g band color, they further argued for low metal- 
licity of its constituent stars. This is consistent with 
the spectroscopic determination by Duffa u et al.l (|2006[ ) 
([Fe/H]=— 1.9), who identified a velocity peak of RR 
Lyrae stars in the same direction, naming it the "Virgo 
Stellar Stream" (VSSfl. Subsequent kinematic studies 
showed that there ma y be individ ual filaments within 
the VOD/ VSS region (!Vivas et all [1)08) . likely to be- 
long to more than one stream. The current knowledge of 
VOD/VSS properties is summarized in Tabled] 

Despite these extensive studies, there is still no defini- 
tive answer on the origin of the VOD/VSS nor of its 
extent. Since the stars are metal poor, but cover a 
wide range o f meta llicities, both Duffau et al. (2006) and 
lJuric et al] (|2008[ ) have argued it to be the debris of 
a tidally disrupted dwarf spheroidal galaxy. S upport- 
ing this hvpothesis. lCasetti-Dinescu et al.l (|2009D showed 
that orbit of an RR Lyrae star associated with the VSS, 
derived from proper motion and radial velocity observa- 
tions, is consistent with a merging scenario in its early 
phase. An alternative, which would explain the unusu- 
ally large angular size of the VOD, is that the overdensity 
may not be due to a merger but a misinterpretation of a 

^ In this paper, we use the term "Virgo overdensity" to denote 
both the excess stellar number density, as well as the kinematic 
peaks observed in the general directon of Virgo at distances ~ 5— 25 
kpc as it is unclear at present whether the two are truly distinct. 



more complex large-scale structure of the Galactic stel- 
lar halo. However, given the limited extent of wide-field 
data available at the time, this hypothesis could not be 
tested conclusively. 

Recently, the Eighth Data Release (DR8) of th e Sloan 
Digit al Sky Survey has become public ( Aiha ra et al.l 
1201 ll ). Compared to previous releases, DR8 adds a sig- 
nificant new area (~ 2000 deg^) in the southern Galac- 
tic hemisphere, as well as filling in a few "holes" in the 
northern footprint. Importantly, the added area allows 
us to conduct symmetry studies, model-free comparisons 
of stellar number densities along directions where we ex- 
pect them to be the same if the halo conforms to a given 
shape (spherical, ellipsoidal, etc.). 

In this paper, we use this newly available dataset to 
remeasure the size of VOD, better understand its large- 
scale structure, and attempt to discern between the two 
proposed options for its origin. We begin this analysis 
in Section [2] with an overview of the used stellar sample, 
details on the photometric parallax, and construction of 
stellar number density maps. Section |3] brings analysis of 
density maps, with special attention given to the extent 
of the overdensity and how it compares to surrounding 
area without overdensities. In the final Section we discuss 
the implications of presented results on the nature of 
Virgo. 

2. DATA AND METHODOLOGY 

In this section we describe characteristics of the SDSS 
imaging survey and the subset of its stellar sample used 
in this work. We discuss how the distance to each star 
was calculated using the photometric parallax method, 
paying special attention to metallicity effects. In conclu- 
sion, we show how the stellar density maps were created. 

2.1. SDSS DR8 Imaging Survey 

The Sloan Digital Sky Survey III provides an 
unabridged view of the night sky. Its total imaging foot- 
print covers 14,555 deg^, a third of the celestial sphere, 
including ^^5,200 deg^ of imaging on the So uthern Galac- 
tic H emisphere available in Data Release 8 (jAihara et al.l 
120 lit ). All of the planned imaging is now complete, but 
spectra will continue to be taken until the project ends 
in 2014. These will further enlarge the existing library 
which already contains spectra of 520,000 stars, 860,000 
galaxies and 120,000 quasars. For a detailed overview 
of ongoing spectroscopic projects see [Eiscnstcin et al.l 

(Hon). 

The SDSS photometric component collects data in five 
optical bands: it, r, i arid z measured in Ai?^, mag- 
nitude system (iGunn et al.l 119981: iFukugita et~all 119961: 
IGunn et al.ll2006l ). with 95% completeness levels at mag- 
nitudes 22.1, 22.4, 22.1, 21.2, and 20.3 respectively. 
The latest data release has very accurate purely inter- 
nal photometric calibration (jPadmanabhan et al.l l2008l 
so-called, iibercalibration), with the same standard-star 
derived zeropoints used to absolutely calibrate the pre- 
vious releases. Furthermore, all the imaging data was 
reprocessed with the new photometric pipelines, featur- 
ing enhanced sky-subtraction algorithm. Unfortunately, 
due to a mistake in data processing, the absolute astrom- 
etry of DR8 is less accurate than the previous releases 
(M. Blanton; private communication), but this is of no 
consequence to the work presented in this paper. 



TABLE 1 

Overview of Virgo Overdensity / Stellar Stream 
properties 



Quantity 


Value 


Refer E^'CE 


Angular size 


> loCI 


iDuffau et a;n (12009) 


(deg^) 








> 1000 


[^uric et al. (200^ 




~ 76(0 


[Prior et al. ( 2003) 




> 2000 


this work 


Surface bright- 


32.5 


Duric et al.l 120081) 


ness (mag 






arcsec ) 






Distance (kpc) 


20 


Ncwberg et al.l 




(2002) 




19 


^ivas & Zinn ( 2003) 




6 to 20 


puric ct al. (2003) 




IS 


prior et al. ( 2003) 




15 to 30 


[Brink et al. ( 2010|) 


Metallicity 


-1.86 ±0.4(0 


Buffau et al. £2006) 


[Fe/H] 








-l.ffl 


lJuric et al.l (120081) 




-2.0 ±0.1 (internal) 


lAn et al.l (120091) 




±0.5 (systematicJ3 




Radial velocity 


99.8 ±17.^ 


IDuttau et al.l (120061) 


(km s"-') 








130 ± 10 


INewberg et al.l 






(120071) 




127 ± 1(0 


[Prior et al." ('2003) 




134.4 ± 14.(3 


[Casetti-Dinescu et al.l 






(f2009i) 


Proper motion 


^aCOsS = —3.50 ± 


ICasetti-Dinescu et al.l 


(mas yr"-') 


0.85, us = 2.33±0.85 


120091) 


Origin 


Sagittarius dSph 


Marti'nez-Deleado et alj 






(|20071) 




dwarf galaxy 


IJuric et al.l (I2008D: 



'Casetti-Dinesc u et al.l 
(2009); this work 



''Quantity related to the VSS. 
''Spectroscopically derived quantity. 
"^Photometrically derived quantity. 

Importantly, the morphological star-galaxy separation 
algorithm is well understood and has not changed since 
DR2. It classifies objects as 'GALAXY' if: 

psfMag - cmodelMag > 0.145. (1) 

This admittedly simple criterion has been shown to work 
well for sel ecting clean samples of stars to r ~ 21.5 
(|Juric et al.ll2008r i. 

2.2. The Photometric Parallax Relation and Iterative 
Determination of Distances 

More than 95% of stars detected by the SDSS are on 
the main sequence and of similar age, therefore resid- 
ing on a fairly constrained, one-dimensional, stellar lo- 
cus. Provided we have a calibrated color-luminosity (or, 
photometric parallax) relation, this fact allows us to esti- 
mate both their absolute magnitudes and distances from 
multi-band photometry. 

Once the absolute magnitude of a star is known, its 
distance is easily calculated using: 

D(pc) = io('--*^'-)/5+i (2) 

which then leads to its position in Galactocentric Carte- 
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sian coordinates: 



X = Rq — D cos I cos b 
Y = —D sin / cos b 
Z = Dsm{b) 



(3) 



where Rq is distance from Sun to the Galactic center and 
/ and b are Galactic longitude and latitude, respectively. 

There have been numerous efforts through studies 
of stellar systems with known distances (eg. nearby 
stars, globular clusters) yielding a number of pro posed 
photometric pa ralla x relations (eg. i Hawlcv ct al. 2002; 
iWihiams et all 120021 : IWest et al.ll2005i : Bilir et al..,2006) 
In this paper, we use the relation derived bv llvezic et al.l 
pOOa : 

Mrig - I, [Fe/H]) = M^,{g ~ i) + AM,([Fe/H]) (4) 

where the terms M^{g — i) and AMr([Fe/H]) have been 
determined to be: 

M°{g ^i) = -0.56 + 14.32(.g - i) - 12.97(.g - if 

+6.127(5 - if - 1.267(.g - if + 0.0967(5 - ifi^) 
A Mr([Fe/H]) = -l.ll([Fc/H]) - 0.18([Fe/H])2 

This relation was calibrated from SDSS observations of 
11 star clusters in metallicity ran ge from +0. 1 2 to ~ 2.50. 
It builds on the previous work of lJuric et all (|2008D . and 
is in good agreement with other proposed relations in the 
literature. It is ex pected to be accurate to 10 — 15% 
(jlvezic et al.ll2008[ ). and the fact that it was directly cali- 
brated on SDSS photometry makes it especially apropri- 
ate for use here. 

While more accurate, the llvezic et al.l (|2008[ ) relation 
is at a disadvan tage compared to, for example, the 
lJuric et aD (|2008[ ) due its need for the u-band photome- 
try to determine metallicity. Given the shallower depth 
of the SDSS u band observations, this would effectively 
limit our explorations to ^ 8 kpc, heliocentric. We there- 
fore chose not to determine metallicity from the pho- 
tometry, bu t instead use the metallicity prior given by 
llvezic et al.1 (pOOS. ) models t o compute the AM rdFe/B]) 
needed in Equation [S] Since llvezic et al.l (|20Q8l ) give the 
metallicity as a function of position in the Galaxy (that, 
in turn, depends on the absolute magnitude), this makes 
our problem an implicit one. 

We iteratively solve for D and Mr starting with an 
initial guess for metallicity of [Fe/H] ~ —0.5. We then 
compute D and Mr from the observed i. Next, the 
expected metallicity at that positio n in the Ga l axy is 
drawn from distributions given by 'Ivezic et al.1 ()2008l 
Equations 18-20), and the process is repe ated until con - 
vergence is obtained. Note that since llvezic et all (pOOl ) 
give their metallicity distributions separately for the disk 
and the halo, we randomly assign a star to the disk or 
halo component, with the weight given by their local con- 
tributions. In practice, this detail is of small importance 
for stars at approximate distance of the VOD, as nearly 
all belong to the halo. 

We have compared the stellar number density depen- 
dence on the position in the Galax y as obtained with 
the iterative approach, to results of llvezic et al.l (I2008D 



as well as the result obtained using lJuric et al 
photometric parallax relation. We found them to be in 
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agreement to within ~10% in the thin and thick disk re- 
gions, and discrepant with respect to the normalization 
of the halo component. In particular , while the overall 
halo density profile given bv lJuric et al.l (|2008l ) is correct, 
we have found their local halo-to-disk normalization to 
be over-estimated by approximately a factor of three. 

For Galactic models used in subsequent sections, we 
have lowered the ha lo-to-thin d i sk no rmalization pa- 
rameter fn from the lJuric et all ()2008D value to fn = 
1.6 X 10"'^. The difference is not sur prising in light of a 
relatively large error bar (50%) that lJuric et alj (|2008f ) 
have attached to their measurement of /h- We also note 
that the value used here is more consistent with those 
traditionally used. A formal fit for the fn is beyond the 
scope of this work and will be discussed in a subsequent 
paper (Juric et al.; in prep). In here, we will refer to it 
as the galfast model. 

2.3. Sample Selection and Maps of Stellar Number 
Density 

To be included in our sample of main sequence stars, 
we require the objects to be morphologically classified as 
'STAR' by the SDSS pipehne, be within ~ q.32 mag of 
the st ellar locus as defined by Equation 4 of lJuric et al.l 
(|2008D . and have colors g — r > 0. The last cut re- 
moves extremely blue point sources, such as BHB stars 
or blue stragglers. We further require the objects to 
have r,i < 21.5, to restrict ourselves to the range where 
SDSS' morpholog ical star-galaxy separation works well 
(jJuric et al.l l2008f l . The apparent magnitudes of the 
resulting sample, containing ^ 86 million objects, are 
then corr ected for interste llar extinction using the SFD 
maps (Sch legel et al.lll998l ). Finally, we apply the iter- 
ative distance/absolute- magnitude determination proce- 
dure to all objects (stars) of this sample. 

We proceed to divide up the sample in 6 absolute mag- 
nitude bins (spanning 3.5 < Mr < 9.5), and 9 shells 
of heliocentric distance (spanning 2 < D/kpc < 20). 
The bin sizes in absolute magnitude and distance are 
1 mag and 2 kpc, respectively. Results shown here con- 
tain only stars from the brightest absolute magnitude bin 
(3.5 < Mr < 4.5), which probe deepest into the stellar 
halo. 

These data are then binned spatially, in Lambert equal 
area projection, and plotted in hemisphere plots. The 
projection poles were set to the north and south Galac- 
tic pole for examination of northern and southern SDSS 
data, respectively. The Lambert projection was chosen 
as it allows for straightforward comparison of sizes of ob- 
served structures on different parts of the sky. The map 
pixel scale of dx = 1.5° was chosen to keep most of the 
pixels well populated, thus reducing Poisson noise. At a 
central-distance bin of 9 kpc ~70 % of pixels contain at 
least ^-^30 stars, while at the largest distance of 19 kpc 
this fraction is '^50 %. 

Not all pixels on the map have been fully covered by 
SDSS observations. This is particularly true for pix- 
els close to the edge of the footprint, and needs to be 
accounted for. For consistent comparison of data and 
model predictions, we compute the pixel fraction covered 
by the SDSS by subdividing each pixel with a 100x100 
grid, and counting the number of so defined subpixels 
that contain at least one star. The fraction of pixel area 
covered by the SDSS is then approximated by the ra- 
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tio of counted subpixels, to the total number of sub- 
pixels. Model predictions (discussed below) have been 
multiplied by this fraction. 

Examples of stellar density maps thus obtained are 
given in Figures [T] (north Galactic hemisphere) and [2] 
(south Galactic hemisphere). The panels in the left col- 
umn show SDSS DR8 stellar counts, while galfast model 
predictions for the same population of stars are given in 
the center. Model predictions give the expected number 
of stars in the volume determined by pixel and distance 
bin size. Finally, a quantitative comparison of the data 
and model, given by their difference normalized to the 
model, is shown on residual plots in the right column. 

3. RESULTS 

In this section we discuss the density maps constructed 
as described above, with emphasis on the structure of 
Virgo Overdensity. Previous studies found the VOD to 
peak at '-~^10 kpc. Here we present and discuss density 
distribution maps at two distances: 5 ± 1 kpc where we 
do not see signature of the VOD, and 11 ± 1 kpc where 
it is very prominent. 

3.1. Analysis of star counts maps 

As seen in Figure [TJ the 5 kpc density map is well 
matched by the model prediction, resulting in mostly 
uniform residual map centered around zero. The density 
map itself reveals the bulk structure of Milky Way and 
our position inside it. The greatest density is observed 
at low Galactic latitudes (6 < 30°) where our distance 
shell dips into the Galactic disk. In order to reduce the 
clutter and obtain a better view of the halo stars, the 
plots show only the region b > 30°. In comparison to 
the model predictions, we note two features on the map 
of residuals. First, the slight underdensity in the region 
heavily populated by disk stars (0° < / < 30°, b < 45°) 
due to model parameters tuned to best reproduce the 
halo, and second the Monoceros stream viewed as an 
overdensity on the opposite side of the map (90° < / < 
230°, b < 40°). 

The situation is radically different at 11 kpc. Even 
though the general features seen in the 5 kpc plots are 
also present at 11 kpc, the residuals are dominated by a 
large overdensity at / ~ 300°, b ~ 75° - the Virgo Over- 
density. Probability of a random fluctuation expanding 
over such an area is statistically insignificant. 

We also present the stellar density maps of the south- 
ern Galactic hemisphere derived from SDSS DR8 data, 
in the same distance shells as for the northern hemi- 
sphere (Figure [2]). The SDSS footprint of southern 
sky, limited by geographical location of the telescope, 
is smaller than in the north, but still large enough to re- 
solve structures on scales as large as several hundreds 
deg^. Apart from the known Hercules- Aquilae cloud 
{I ~ 60°, b ~ —35°) and trailing arm of the Sagittarius 
dwarf galaxy (90° < I < 230°, b > -40°), there appear 
to be no strong overdensities in the halo region. How- 
ever, we do observe a curious density enhancement in the 
south at low latitudes in the anticenter direction. This 
overdensity, that could be associated with the Monoceros 
stream, or be a feature of the thick disk, will be analyzed 
in a subsequent publication. Other than the mentioned 
overdensities, the residual maps exhibit a very good fit 
to the model. 



3.2. Extent of the Virgo Overdensity 

The fact that the VOD is not only visible on the resid- 
ual plot, but also as an enhancement on the density plot 
itself, motivates us for another, model free look at the 
data. This approach also reduces the dependence of 
our results and conclusions on the details of the analytic 
Galactic model. 

Most overdense pixels on the VOD residual plot (lower 
right panel on Figure [Ij are above the I = 0° line. If 
the halo were symmetric, the halves of the plot divided 
by this line should also be symmetric. In Figure [3] we 
compare the stellar density in these halves by plotting 
the value of: 



fw-E{l,b) 



data{l, b) 



data{360° -l,b) 



1 



(6) 



For a given pixel at Galactic coordinates {I, b), the value 
of fw-E is simply the normalized difference between the 
star counts in the pixel itself (data{l,b)) and its coun- 
terpart (data(360° — l^b)). The quantity fw-E is de- 
fined only for the pixels with axisymmetrical counter- 
parts, which results in the "bug like" shape of the plots. 

The stellar number density was found to be symmetric 
on the order of < 10% at distances closer than 5 kpc. 
However, at 11 kpc the median of fw-E distribution 
on the upper half is at ^ 20%. If we adopt the def- 
inition of overdensity as density contrast (here defined 
with respect to the axisymmetric region) of >15%, then 
the total area of Virgo overdensity as measured in SDSS 
pR8 is » 2000 de g^, double the size previously measured 
(jJuric et al.l[2008l ). 

Furthermore, we see that the overdensity extends along 
the I = 300° and 270° SEGUE stripes. This indicates 
that 2000 deg^ is still only a lower limit of its actual 
size. We attempt to roughly estimate the area of the 
VOD outside the SDSS footprint by placing an ellipse 
encompassing most of the pixels with 40% overdensity in 
the SEGUE stripes as well as the overdensity in the main 
area of the survey. We find this area to be ~3000 deg^ , 
implying that approximately one third of the VOD may 
extend over the area where no data is currently available. 

3.3. Luminosity of the Virgo Overdensity 

To further characterize the progenitor, we estimate a 
lower limit on the number of stars associated with VOD. 
The Virgo area was conservatively limited to region sub- 
tended by: (200° < I < 210°, b > 40°), (210° < / < 270°, 
5 > 50°) and (270° < I < 340°, 6 > 60°). To estimate the 
excess in stellar number counts, we simply subtract the 
pixel values in symmetric, Virgo-free area from the VOD 
pixels, and sum over all the pixels and distance range (7 
- 19 kpc). This gives the number of VOD stars w ith ab- 
solute magnitude Mr — 4.0 ± 0.5. Following Juri c et ahl 
(2008:), and assuming the luminosity function of VOD 
stars is similar to that of the halo (Juric et al., in prep), 
we obtain an order-of-magnitude estimate of ^ 10^ for 
the total number of stars present in VOD. This is consis- 
tent with numbers found for large globular clusters and 
dwarf galaxies. 

3.4. Test of the "Tilted Halo" model 

Given the large extent of the VOD, an attractive hy- 
pothesis is that it may not be an overdensity in the stellar 
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Data 



Model 



Residual 




Fig. 1. — Density distribution of F stars in 2 kpc wide shells centered on 5 and 11 kpc, mapped in Lambert equal area projection with 
circles representing constant Galactic latitude and lines constant Galactic longitude. Disk area is excluded by plotting only latitudes higher 
than 30°. Left p anels show SDSS DR8 imaging data, with blank areas corresponding to pixels with no data, best-fitting models updated 
from lJuric et al'l (2008) for the same area, corrected for SDSS sky coverage, arc on the middle panels, while the data - model residuals 
(normalized to the model) are on the right panels. At the 5 kpc distance shell, data and model are in very good agreement (apart from 
the disk region, which was not modeled as carefully), but the 11 kpc shell clearly shows Virgo overdensity extending over 2000 deg"^ in 
general direction of I ~ 300° , b ~ 75° . The Monoceros stream is visible at both distance shells as overdensity on low Galactic latitudes for 
150° < I < 230°. 



halo but a signature of a more complex stellar halo den- 
sity distribution. For example, proposals have been put 
forward that the Hercules-Aquilae cloud is a signature 
of dynamical interaction of the disk with the stellar bar, 
and not a merger event (Hum phreys et al.ll20ll[ ). 

An oblate halo model, with axes aligned with those 
of the Galactic disk, does not produce an overdensity 
signature in maps constructed in Section 2.3. However, 
if its axes were not aligned with the Galactic disk (eg. 
if it is "tilted"), such an overdensity signature will 
occur (Figure [J) . With the SDSS data in the southern 
hemisphere available for a significant portion of the sky, 
it is possible to test this tilted halo model by directly 
comparing the stellar number density above and below 
the Galactic plane. 

We proceed similarly to steps taken in previous section 
and plot the values of /n-s defined by: 



fN-s{l,b) 



data{l, b) 



(7) 



data{l, —b) 

in Figure [5l data{l, b) once again denotes stellar number 



count in pixel centered at (/,6) and distance of 11 kpc, 
analagous to the bottom rows of Figures [1] and [21 Note 
how the value of fN-s is only defined for those pixels 
which have data available both above (6 > 0) and below 
(5 < 0) the Galactic plane. Although the distribution of 
fN-s values is noisy, there is no clear over- or underden- 
sity as would be predicted by Figure IH This is further 
confirmed by the distribution of densities in pixels not 
affected by overdensities (excluded regions include Sagit- 
tarius dwarf galaxy at |6| < 60° and 120° < I < 150°, 
45° < |6| < 60° and Hercules-Aquilae at circle of ra- 
dius 7° around I = 55°, 6 = -30°). The distribution is 
slightly asymmetric, but the median is at -0.01, while the 
width of the distribution (with semi-interquartile range 
of 0.22) is consistent with the noise visible on the map 
of /at-s values. Given the distribution centered around 
zero excess of north versus south, we can conclude that 
on the scales of ~20 kpc, halo is symmetric with respect 
to the Galactic plane. 

This remains true when the exercise is repeated for dis- 
tance shells in 7 — 19 kpc range. Figure [6] shows how the 
median of /n-s distributions similar to that on Figure [5] 
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Fig. 2. — Same as Figure[T]on southern Galactic hemisphere. Apart from the enhancement at low latitudes in the anticenter region that 
may be related to the Monoceros stream, the Sagittarius stream (90° < I < 230°, b > —40°), and the Hercules- Aquilae cloud (/ ~ 60°, 
b ~ —35°), no new strong excesses were detected. 



changes with distance. The error bars denote the distri- 
butions' IQR. The plot shows remarkable consistency in 
difference of star counts between northern and southern 
Galactic hemisphere, the largest deviation being ^ 10% 
at 7 ± 1 kpc, and considerably smaller in most other dis- 
tance bins. 

We therefore conclude that a tilted halo cannot ac- 
count for the existence and observed extent of the Virgo 
Overdensity. 

4. DISCUSSION 

Previous studies determined the distance range, radial 
velocities and metallicity of stars making up the VOD 
(see Table [Ij . Several authors have also reported evi- 
dence of substructure (e.g. Duffau ct al. 2006). In this 
paper we have used the SDSS DR8 data to reasses the 
extent of the VOD and find it likely to extend over 
^ 3 000 deg^ of the sk y, three times more than the origi- 
nal [JuriHeril] (|2008f ) estimate. 

Stellar overdensities are usually attributed to merger 
remnants, but VOD's size and morphology make it un- 
like a typical tidal stream thus calling for a more cau- 
tious interpretation. An alternative to the merger rem- 
nant hypothesis is that the VOD is a signature of a more 
complex global halo density distribution (for a discus- 
sion why oth er pro posed hypotheses are less likely, see 
section 5.5 in lJuric et al...2008. ). In this paper, we have 



tested whether a "tilted" halo can explain the observed 
signature. 

Figure S] shows schematically how a tilted halo would 
appear as an overdensity in Figures [3] and [51 Axisym- 
metric halo on the left panel would be observed as being 
completely symmetric in both fw-E and Jn-s values 
defined earlier, but if the halo were tilted as seen in the 
right panel on Figure IH a VOD-like enhancement in stel- 
lar number counts would be observed in the "western" 
side of the fw-E plot. Halo tilt would also be evident on 
the Jn-s plot because the number counts on the "south" 
would be larger compared to the symmetric area on the 
"north" . However, as we can see from Figure |6l this is 
not the case. In fact, given the difference in stellar num- 
ber counts between north and south areas mapped by 
the SDSS, we can constrain any north-south halo asym- 
metry to be less than 10%. As the density contrast of 
the VOD gets as high as ^^50% (see Figure [5]), it is clear 
this simple tilted halo model can only amount to a minor 
contribution to the overall enhancement in density. This 
leaves the merger hypothesis as the most likelj^. 

While the VOD's morphology is clearly different from 
those of other streams (eg. Sagittarius, GDI, etc), 

^ While one can always invoke even more complex halo shapes to 
explain the observations, at some point these become unreasonably 
contrived. 
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Fig. 3. — A model free confirmation of the Virgo Overdensity 
is given by a star count difference between 0° < I < 180° and 
180° < I < 360° regions on the bottom left panel from Figure [l] 
The above region is ~ 50% more densely populated than the lower 
region. The lower limit on the overdensity area is ~2000 deg"^, 
but note that the overdensity extends down to 6 ~ 30° along the 
I = 270° SEGUE stripe, so the real extent of the Virgo Overdensity 
might be considerably larger. 



N N 




Fig. 4. — Schematic view of the cross section of oblate halo (dark 
gray) and spherical shell centered on the Sun (light gray). The 
axis connecting the Sun and the Galactic center is perpendicular 
to the surface of the plot. In the case of axially symmetric halo (left 
panel), the stellar density counts at certain distance (as marked by 
the light gray shell) are symmetric with respect to both I = 0° line 
(W vs. E) and the Galactic plane (N-S). However, if the halo is 
tilted (right panel), so that we observe excess on the "Western" 
part of the sky (ie the Virgo Overdensity), then the North - South 
symmetry holds no longer either. 

its cloudlike appearance is not entirely unique. Other 
such overdensi ties have been observ ed, eg. Hercules- 
Aquila cloud (IBeloku rov et al.l 120071) . or Triangulum- 
Andromeda plocha-Pi nto et all I2004D . Similar struc- 
tures appear in Gala ctic halo fo rmation simulations as 
well. For example, the lJohnston et al.. (,2008. ) simulations 



Aihara, H., et al. 2011, ApJS, 193, 29 
An, D., et al. 2009, ApJ, 707, L64 
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— . 2007, ApJ, 657, L89 



modeled to match observed properties of Local Group 
galaxies predict several cloud-like morphologies with sur- 
face brightness Ri32.5 mag arcsec"^ in a Milky Way type 
galaxy. This is broadly consistent with what we see in 
Virgo Overdensity; however, we do find traces of VOD 
as close as 7 kpc from the Galactic center, closer than 
the simulations would predict. Given the simplicity of 
the simulations, this may not be a serious issue. 

Furthermore, the accretion origin of the VOD is con- 
sistent with kinematic measurements of an RR Lyrae 
star assumed to belong to the Virgo Stellar Stream 
(jCasetti-Dinescu et aLH 2009). Based on the radial veloc- 
ity and p roper motion measurement w ith a baseline of 
a century, iCasetti-Dinescu et al.l ()2009D have determined 
the orbit of an VSS RR Lyrae, having Vperi = 11 ± 1 
kpc, rapo = 89tg kpc, and a period of T = 1.21°;^ Gyr. 
The measured orbit passes through the VOD, indicating 
that VOD and VSS may be connected structures with a 
common origin. Also, the high eccentricity of the orbit 
is consistent with the observed high eccentricities in sim- 
ulated merger events that result in cloudlike structures. 

Taken together, all these lines of evidence point to 
the VOD being a signature of a high eccentricity merger 
event observed at perigalacticon. Future kinematic 
studies are the best way to further verify this conclusion. 
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Fig. 5. — A map of difference in star counts in nortfiern and soutfiern Galactic hemisphere (see bottom left panels on Figures [T1 and [2l 
respectively), normalized to the south. A very good match of the number counts over 20 kpc range provides strong constraints on the halo 
symmetry with respect to the Galactic plane and disfavors the tilted halo interpretation for the Virgo Overdensity. 



0.4 - 



0.2 - 



"O 

-0.2 - 



-0.4 - 



8 10 12 14 16 18 

Distance (kpc) 
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Equation O is here computed only for the areas not contaminated by Hercules-Aquilae cloud and Sagittarius stream. The error bars 
represent the semi-interquartile ranges of the /jv-S distribution. The biggest deviation from symmetry (denoted by the dashed line at 
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